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Abstract 

The canonical quantum theory of gravity - Quantum Geometro¬ 
dynamics (QG) is applied to the homogeneous Bianchi type IX cos¬ 
mological model. As a result, the framework for the quantum theory 
of homogeneous cosmologies is developed. We show that the theory 
is internally consistent, and prove that it possesses the correct classi¬ 
cal limit (the theory of general relativity). To emphasize the special 
role that the constraints play in this new theory we, compare it to 
the traditional ADM square-root and Wheeler-DeWitt quantization 
schemes. We show that, unlike the traditional approaches, QG leads 
to a well-defined Schrodinger equation for the wave-function of the 
universe that is inherently coupled to the expectation value of the 
constraint equations. This coupling to the constraints is responsible 
for the appearance of a coherent spacetime picture. Thus, the physical 
meaning of the constraints of the theory is quite different from Dirac’s 
interpretation. In light of this distinctive feature of the theory, we 
readdress the question of the dark energy effects in the Bianchi IX 
cosmological model for highly non-classical quantum states. We show 
that, at least for this model, for any choice of the initial wave function, 
the quantum corrections will not produce the accelerated expansion 
of the universe. 
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1 Introduction 


The canonical method of quantization is one of the most direct paths to the 
quantum theory. It is appealing because it is based on the well-established 
classical theory, and therefore leads to a plausible physical theory by the 
virtue of Erenfest’s theorem. Therefore, the resulting quantum theory auto¬ 
matically reproduces, at some limit, the observable physical effects described 
by the classical theory. This feature is true for the general quantum gravita¬ 
tional field. 

It has been a long-standing problem to apply this method to Einstein’s 
geometric theory of gravity - the theory of general relativity. The challenge 
is to find an appropriate time-translation operator, a Hamiltonian of the the¬ 
ory, which would govern the dynamics of the theory at the quantum level. 
The problem with general relativity is that the Hamiltonian is the genera¬ 
tor of gauge transformations and therefore must vanish. Simply speaking, 
it is a constraint of the theory. This does not cause any problems for the 
classical dynamics, since one can still calculate the Poisson brackets using 
this Hamiltonian and obtain Einstein’s equations of motion. However, at the 
quantum level, according to Dirac’s treatment of the constraints, the operator 
corresponding to the classical expression for the constraint must annihilate 
physical states of the theory. This means that classical Hamiltonian when 
expressed as a quantum operator fails to generate time translations. A chal¬ 
lenge facing canonical quantum gravity is to find an alternative operator for 
time translation, this is often referred to as “problem of time. ” HEain. 

In this article we apply the formalism of Quantum Geometrodynamics 
(QG) 0 El El E| to the anisotropic homogeneous Bianchi type IX cosmolog¬ 
ical model. Because of its structure, it can be parameterized with only few 
degrees of freedom instead of the infinity of degrees of freedom of the field 
in the general case. On the other hand, due to its anisotropy, the model is 
still rich enough to reflect some of the characteristic features of the cosmo¬ 
logical singularity, m and thus may be an interesting candidate to test the 
quantum theory. 

We have two goals in applying the general formalism to this concrete 
system. First, we would like to emphasize the special role played by the 
constraints in this theory, which is conceptually and interpretationally dif¬ 
ferent from the more traditional approaches. Second, we wish to develop 
a quantum cosmology theory, within which, one can address the issues rel¬ 
evant for cosmology in general. In Sec. 2 of the article, we briefly review 
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the general theory ( for extended discussion see |7J) mostly for the sake of 
introducing notations. In Sec. 3 we introduce the Bianchi IX cosmological 
model as a homogeneous solution of Einstein’s equations of motion. In the 
next two sections, Sec. 4 and 5, we apply the ADM square-root quantization 
and Wheeler-DeWitt theory to the model, with the aim to illustrate the rel¬ 
ative successes and failures of these formalisms. We place special emphasis 
on the treatment of the constraints in these theories. In Sec. 6, we develop 
the quantum theory of the Bianchi IX cosmological model within the general 
formalism of QG. Here we show that, following the prescription, it is possible 
to define a Hamiltonian operator which leads to a Schrodinger equation cou¬ 
pled to the constraint equations. We show, quite generally, that the resulting 
quantum theory reproduces Einstein’s equations in its classical limit. The 
detailed proof of this for any system with finite degrees of freedom is given 
in the Appendix A. Within this model we show that the quantum effects do 
not lead to an acceleration of the expansion of the universe. Finally, in the 
Appendix B, we collect the technical details of the invariant basis for Bianchi 
IX spacetimes. 


2 Quantum Geometrodynamics — the Gen¬ 
eral Theory 

We start with the Hilbert-Einstein action for general relativity in its ADM 
or 3+1 formulation in. As usual we choose a foliation of the 4-dimensional 
manifold by a one parameter family of the spacclike hypersurfaces, X t . We 
assume that the spacetime we are solving for is globally hyperbolic, and that 
this kind of foliation can always be accomplished. The canonical degrees of 
freedom are taken to be the six components of the spatial metric hij induced 
on the hypersurface Ik by the full spacetime metric g /w restricted to act on 
vectors tangent to the spatial slice. The spacetime metric tensor then can be 
written as 

Q^iv hfiv T +/ +v: (l) 

where h M are the components of the timelike unit vector normal to the space- 
like hypersurface, X t , in some coordinates {x a }. In what follows, the partic¬ 
ular choice of coordinates we use are irrelevant. However, in order to simplify 
our calculations, we assume that, at every point of the 4-dimensional mani¬ 
fold, the coordinates {x*}, i = 1,2,3, define the vectors {W} tangent to Y, t . 
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The fourth basis vector has to be timelike but otherwise arbitrary, so that 


= Nn M 


( 2 ) 


where N is an arbitrary function on the manifold, and A+ is the vector 
tangent to In the {X i,h} basis, its components ( N,N l ) are the familiar 
lapse and shift functions. This vector defines the flow of time, and since 
we leave N and N l unspecified, it manifests the freedom of choice for the 
timelike direction. With the foliation £ t and basis of tangent space given, 
we cast the Hilbert-Einstein action in its 3+1 form. 


S=^—J VhNf 3 >R+K ij K‘ i -K 2 ]d 3 xdt= +- 

(3) 

Here, K l0 is extrinsic curvature and is defined as the Lie derivative of the 
3-dimensional metric with respect of the unit normal vector n to K is the 
trace of Kij, and h is the determinant of h^. The Legendre transformation 
is accomplished with respect to the “time” derivative of the spacclike metric 
and is defined as, h^ = h i k hj k C t hki. The canonical momenta are then found 
to be 

7 T ij = = Vh(K ij - Kh lj ), (4) 

Shij 

C° = \/h(- (3) /f — h^ir^nij — -/i _1 7r 2 ), and (5) 

C j = -2 VhDi(h- 1/2 TT ij ). (6) 


By varying © with respect to N and N l , we get the constraints of the theory 
which are just C a . By varying the action with respect to the other canonical 
variables and their conjugate momenta, we arrive at Hamilton equation of 
motions, 


hij 


SHg 


2h- 1/2 N(7iij - 


1 

2 


hij 7r) + 2 D(iNj), 


7T 


i] 


SH g 

5 hij 


(7) 

( 8 ) 


Where Hq = J(NC° + NiC l )d 3 x is Hamiltonian of general relativity. 

Equation m is just a definition of the momenta expressed in terms of 
the “time” derivative of the metric of the 3-dimensional hypersurface. It is 
a good place to comment on the latter. In most of the sets of basis vectors 
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used in practice, this definition of the time derivative, which is based on the 
concept of the Lie derivative, becomes just a partial derivative with respect 
to the coordinate time, t. The conditions for this is that t commutes with 
the spacelike vectors, X *, of the basis. In other words it must be coordinated 
with respect to the spatial basis. This condition is trivially satisfied in the 
case of the basis defined by the coordinate functions on the manifold, 

but one has to make sure that it still holds in case one decides to work in 
some other basis, say invariant basis as for Bianchi IX model we consider 
here. 

When we have the classical theory cast in its canonical form, it can be 
quantized. The key idea of the new formalism is to split the classical canoni¬ 
cal degrees of freedom, which are the six components of hij, into the so called 
true dynamical variables ({<&}, d — 1 — 2) and the embedding variables (g e , 
e = 1 — 4). Moreover, we treat only the true dynamical variables {g^} as the 
genuine physical degrees of freedom of the gravitational field, and as such 
are susceptible to quantum fluctuations. Whereas, the set of four embedding 
variables are thought of as classical parameters of the theory, and responsi¬ 
ble for the appearance of the coherent spacetime picture. Following York’s 
analysis of canonical structure of general relativity, we identify the physical 
degrees of freedom with the conformal 3-geometry of the spatial slice. m 
Next, we define the dynamical Hamiltonian of the theory H dyn (p d ,q d ,q e ,q e ) 
by performing the Legendre transformation with respect to the true dynam¬ 
ical variables only. Consequently, the standard canonical methods of quan¬ 
tization lead us to a Schrodinger equation for the functional T( q d ), where in 
the expression for H dyn , q d and p d becomes linear operators with the standard 
canonical commutation relations. 



(9) 


To make the system mathematically complete, we specify the values of the 
embedding variables on each S t . This is achieved by imposing the constraints 
of the theory as expectation value equations. 


(C l (p d ,q d ,q e ,qe)) = 0 
(C°(p d , q d , g e , Qe)) = 0 


( 10 ) 

( 11 ) 


This gives us a set of four first-order differential equations for the four 
unknown embedding variables. Note that these equations do not impose 
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restrictions on the Hilbert space of the system, but instead can be thought 
of as an implicit definition of the embedding parameters. 

Therefore, in QG the time-dependent Schrodinger equation for the state 
of the gravitational held is coupled to the set of expectation-value constraint 
equations. This feature of the theory seems to be as intriguing as it is un¬ 
avoidable, and must result in observable effects when applied to a realistic 
problem. 


3 Bianchi IX Cosmology 

In this section, we introduce the notation and conventions we use to describe 
the homogeneous Bianchi IX cosmological model. 

The characteristic feature of the Bianchi IX model is the existence of the 
simply transitive isometry group jT5|. The infinitesimal generators of this 
group are the three linearly independent spacclikc killing vectors, which 
obey 

[&.&]= CV*’ (12) 

where structure constants are C k l3 = £kij- 

The typical orbit for this group of a point of the four-dimensional man¬ 
ifold, is a spacclikc hypersurface. One can fill the entire space with a one- 
parameter family of such hypersurfaces (£ t ), so that the manifold can be 
represented as a direct product of the real line and this family of hypersur¬ 
faces, St x R. We use the parameter which labels hypersurfaces as the fourth 
coordinate together with other three that covers the spacelike sections. For 
the latter it is convenient to use the coordinates of the three-sphere 0 < 9 < 7r, 
0 < (p < 2-7T, 0 < if; < 47T due to the fact that the isometry group of this 
model is isomorphic to the 3-dimensional rotation group. 

Given the Killing vectors {£*} in these coordinates, one can generate the 
so-called invariant basis {t, e*}. The details of this procedure can be found 
in Appendix B. In this new basis the metric is 

g = —(N 2 - N r N r )a° 2 + N i( r 0 a j + %(£)c tV, (13) 

where N t = hijN and a a is the dual basis of the {£, e*}. 

The advantage of the invariant basis is that we are using the symme¬ 
tries of the spacetime to get rid of the explicit dependence of the metric on 
spatial coordinates. In what follows, we work in this more general invariant 
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basis in order to keep the lapse and shift functions (N(t) and N l (t )) unspec¬ 
ified. That allows us to mimic the analogous degrees of freedom in the full 
Hilbert-Einstein action without special symmetries. On the other hand, we 
must point out that we do use the particular slicing of the four dimensional 
manifold with the assigned Bianchi type IX isometry group. This fact might 
not be considered a great disadvantage from the “3+1“ point of view. This 
slicing provides a natural way to build the spacetime possessing the required 
symmetries as a time evolution of a 3-dimensional space. 

The metric on the spatial slices, X 4 , of the vacuum Bianchi IX model can 
be parameterized in the following manner: 

Kj = R 2 0 e 2n (e 2/3 ) _ a 1 a 3 ; (14) 

where Rq is initial radius of the universe, and all other parameters are func¬ 
tions of time only. The matrix 

flij = diag[/3 + + \/3/3_, /3+ - V3/3-, -2 f3 + \, (15) 

with the property Trj3 = 0, ensures that the 3-volume of the hypersurface 
depends only on the conformal factor 11. In particular the 3-volume of the 
universe is given by V universe ~ R$e 3Q . 

With the use of the metric, one can calculate the 3-dimensional quantities 
necessary for building the action. 

(3, ' R = ^( 1 -' / (/3 + ,/3-)) (16) 

ZKo 


V(P+, (5-) = i Tr (e 4/3 - 2e“ 2/3 + l) = ^e 4/3+ (cosh 4^3/3- - 1) + 1 

— ^e~ 20+ cosh 2v / 3/?_ + \e~ 8p+ (17) 

o o 

Vh = R 3 0 e 3n (18) 

Taking into consideration that all the canonical variables, including the 
lapse and shift functions, are independent of the spatial coordinates, and 
that this model satisfies the momentum constraints identically, we get 


S = n 


dt 


f)h • ■ 

~' 3 -^ - NC° 
dt 


(19) 
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Here we have used the fact that / a 1 A a 2 A a 3 = (47r) 2 . By using expressions 
for the momenta we derive the first Hamilton equation. 


dhij 

~df 


_ K ^ x 
2 V 17 k ' ’ 


( 20 ) 


This can be solved for the momenta in terms of time derivative of the metric. 


7r. 


y/h . 


13 ~ 2 N^ 3 


hijh\) 


( 21 ) 


Here, h k k = h lj hij. 

The special slicing and convenient choice of basis we adopt has allowed 
us to reduce our problem to the dynamics of a system with only a finite 
number degrees of freedom. We further simplify the model by substituting 
the particular form of metric we have chosen in to the general expressions. 


dhij 

~dt 


20 hy T hij fiij 


(no summation!) 



7T 


ij dhij 

dt 


= 12 ; ^(-S7 2 + 0+ 2 + pj) 


C° = (- (3> fl+ 2)j(-f! 2 + P+ 2 + /3- 2 )) 

Putting the two last terms together we get, 

S = 7! j ^h2[-6f2 2 + 6/3+ 2 + 6/1 12 + N 2{3, R]dt. 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


In principle, this action (ESI) might be the starting point of the theory for 
the Bianchi IX cosmology, but then the relation to the full theory would be 
obscured. In the next three sections we apply (1) the ADM square root, (2) 
the Wheeler-DeWitt and, (3) the QG quantization schemes to the system 
described by this action, respectively. 







4 ADM Square Root Quantization 

The main idea of the ADM quantization scheme is to solve the constraints 
before quantizing, and to use four out of the six canonical variables as coor¬ 
dinate labels. The action for the theory, which must be quantized, become 

Sad Ad = yy— J d 4 x (27) 

The ADM action for the Bianchi IX cosmological model is 

/ • • 774 • 

P+d(3 + + p-d(3- - —QdQ. (28) 

In the expression above we got rid of the explicit dependence on the coor¬ 
dinate time. Following the approach of Misner, |13| one treats Q as the time 
parameter for the evolution of the conformal 3-geometry, then by definition, 
the term in the action in front of dfl is the ADM Hamiltonian for the system. 
According to this prescription, one solves the constraint equation to find the 
expression for the Hamiltonian in terms of the canonical variables and their 
conjugate momenta. Recall that in our case, the constraint equation might 
be rewritten as follows: 

( 29 ) 

hence, 

h adm = P 2 ++P 2 -- ^ i3) R- (30) 

All we need now is to take the square root of this expression to obtain the 
generator for time evolution. 

This last step is not problematic if one analyzes the classical dynam¬ 
ics, because the magnitude of Hamiltonian ( Hadm = is always a real 
number by definition. Then the constraint equation puts restrictions on the 
initial momenta. 

The quantum mechanical system is different. Here we must define the 
operator, Hadm■ This can be done if we find all of the eigenvalues, E n , and 
corresponding eigenvectors, | (p n ), of H\ dm . If all the eignevalues are positive 
( E n > 0), then Hadm = X) y /r E^\p n )(ip n \. In our case, however, not all of 
the eigenvalues are positive because the expression for H\ DM is not positive 
definite. To see that this is the case, let us set up a Gaussian state around 
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/3_|_ = 0, /?_ = 0 with zero initial momenta and calculate the expectation value 
of the H\ dm at Q = 0. It is easy to see that R) = which means that 

(H\dm) = negative. Obviously it is impossible to get a negative 

expectation value of the operator which possesses only positive eigenvalues. 
This shows that Hamiltonian operator is not Hermitian on its full domain 
and therefore the dynamics generated by this operator will not be unitary. 


5 Wheeler-DeWitt Equation 

In contrast to the ADM procedure, PI the canonical theory of quantum 
gravity suggested by DeWitt Pi does not try to reduce the set of canon¬ 
ical variables by solving constraint equations. Instead it treats all of the 
components of the spatial metric on equal footing by promoting them and 
their conjugate momenta to the linear operators acting on the Hilbert space 
of physical states. The latter is defined by the requirement that the physi¬ 
cally admissible states have to be the eigenvectors of the constraint operators 
with zero eigenvalues. So the constraints are realized as restrictions on the 
Hilbert space of the states of the theory. This approach was first introduced 
by P.A.M. Dirac as a general method of canonical quantization of theories 
with first-class constraints. m Einstein’s theory of general relativity falls in 
this category, so it is logical to apply the general prescription of Dirac to it. 
The first non-trivial step in the procedure is to determine the physical states 
of the theory. In the case of gravity, they are the solutions of 

C^phys = 0, and (31) 

C°*phy. = 0. (32) 

The first condition ensures that the wave-functional, v I'(/ij J ), depends only on 
the 3-geometry of the spatial hypersurface. It can be realized by choosing the 
appropriate coordinate-independent representation for the wave-functional. 
For example, one can require T to be a function of the powers of scalar 3- 
curvature (^ R ). Whichever way one chooses to implement this condition, it 
is not expected to give us any dynamics, simply because it depends on the 
state of the field on the spatial slice and gives one no information on what 
the physical state must be off of the slice. It does not map the physical states 
on different T, t hypersurfaces. 
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The idea of DeWitt was that the last constraint equation (|32p. which is 
quadratic in canonical momenta, ought to serve this purpose. One assumes 
that this equation, by itself, contains all the dynamics of the gravitational 
field. The classical expression for this constraint can be rewritten in the 
following way: 

G ijk iTT ij TT kl - Vh®R = 0; (33) 

where G ljk i = \h~ x l 2 ( h ik hji + h,ih ]k — hijh k i ) is the DeWitt’s supermetric. 
It is the metric on the space of 3-metrics which is a six dimensional manifold 
with hij being six independent components of a tangent vector. Note that 
this notation differs from the standard one by the use of the lower indices 
instead of upper ones. In all other respects one can fully rely on the differen¬ 
tial geometry methods applied to the six dimensional manifold and find the 
inverse supermetric, 


G ijkl = [h ik h jl + h il h jk - 2h ij h kl ) , (34) 

which can be used to contract two vectors. An interesting feature of this 
construction is that Giju has signature [—b + + ++], which means that it 
induces a time-like direction on the space of 3-metrics. This fact supports De- 
Witt’s expectation that the Hamiltonian constraint possesses the dynamical 
content of the theory. 

In the quantum theory, with the particular coordinate representation we 
have chosen, the canonical momenta operators become functional differential 
operators, 

ij . <5 
77 = 

0 tlij 

The constraint equation is a Klein-Gordon-type equation, 



{GijkiTT~T7 -^ Vh^Fl\ — 0, (36) 

V ohij oh k i ) 

where the first term due to the signature of the metric, and is a kind of 
Laplace-Beltrami operator in a six-dimensional Riemannian manifold, with 
Gijki being its contravariant metric. Thus we have Wheeler-DeWitt equation 

m- 

From the theory of Klcin-Gordon equation, it is known that the space 
of solutions can be assigned with a naturally conserved inner product. The 
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inner product for Wheeler-DeWitt equation takes the form, 

(*a, *b) = Z [ n -i(KGw^r - ^rG ijkl * b ) dr-1. ( 37 ) 

X V ° h kl Shu ) 

Here Z is a proper normalization factor, the integral is taken in the space of 
metrics over a product of a set of five-dimensional hypersurfaces X(x) defined 
for each point of three-dimensional spatial slice, and dT,^ is the directed 
volume element for each five-dimensional hypersurface respectively. This 
inner product is explicitly independent of the choice one makes for the X(x) 
at each point of the spatial slice. 

These general expressions, 4HD and (EH), of the Wheeler-DeWitt theory, 
when applied to the particular case of homogeneous cosmologies, are sim¬ 
plified tremendously. The main reduction comes from the independence of 
the 3-metric on the spatial coordinates. Thus, instead of the infinite product 
over all points of 3-dimensional spatial slice, we have just one integration in 
the space of metrics for all points. In addition, the functional differential 
operators become simple partial derivatives. 

I 11 case of the Bianchi IX metric, the 3-metric hij is diagonal and there¬ 
fore is determined by only three independent parameters. In the particular 
parameterization of the metric we use, the Hamiltonian constraint is, 

c0 = s(- p h7+i’--T ,I) «) = 0 ’ OS) 

where pn = — is the canonical momenta conjugate to D. Note that 
we treat all of the variables on the same footing, and therefore, we have to 
introduce the momenta for conformal factor. In ridding ourselves of overall 
factor we obtain the Wheeler-DeWitt equation. 

(-&> + m + w - + ^ )R )* M ^- )=0 m 

The coordinate basis {gfy > gf - } is orthogonal, although it is not an 
orthonormal basis with respect to the supermetric. If we choose the hy¬ 
persurface in the space of the metric as the surface of constant D, then 
up to irrelevant normalization factor, the inner product for the solutions of 
Wheeler-DeWitt equation become, 

r ( <9\Eu \ 

= (40) 
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At this point one just needs to solve the resulting equation, and by specifying 
the appropriate boundary conditions, choose among all possible solutions 
those which represents the physical situation in question. 

It is this last step of the procedure that gives rise to a major flaw of 
the theory. In particular, the natural product m is not positive definite 
for a generic solution of the Wheeler-DeWitt equation. Therefore, if one is 
to give probabilistic interpretation to wave-function, T(f2, /3 + , /?_), one will 
have to restrict the space of the solutions to those that have a positive norm, 
they would be positive frequency solutions with respect some natural notion 
of time. For a general spacetime, it is ordinarily impossible to accomplish 
such a subdivision of solutions into positive and negative frequency solutions 
due to the lack of symmetries. One is faced with the dilemma to either (1) 
assign a different inner product to the wave-function or even turn it into 
a field and accomplish third quantization, or (2) to restrict ones attention 
to a particular class of physically relevant spacetimes whose symmetries or 
asymptotic behavior allow such a division and impose a boundary condition 
which generates it. Despite many efforts, neither of the two paths has led to 
a consistent and general quantum theory of gravity. 


6 Dynamical Hamiltonian and the Schrodinger 
Equation 


Having the action for the system at our disposal (12(1 . we proceed with the 
standard canonical procedure to obtain Hamilton’s equations of motion. We 
deviate here from the textbook prescription and perform a Legendre trans¬ 
formation of the Lagrangian only with respect to the variables j3 + and /?_. 
This promotes the two anisotropy variables as the parameters that represent 
the true dynamical degrees of freedom of the spacetime. We determine the 
momenta conjugate to j3 + and (3 _. 


P+ = 


V- = 


dL 


G 


d(3+ 

dLg 

d/T 


127 ri?oe 3Q ■ m ■ 
N ^ = 1W + 


12nRle 3fi ■ m ■ 
N ; ” = j\T~ 


(41) 

(42) 
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We define the dynamical Hamiltonian as 


Ho m = pj + +P J_- La = N^ + ^) + ^ . (43) 

The action, written in terms of the new variables, is thus 

S = f (p+P+ + P-P- ~ H Dyn (p + ,p-,p+, /3_, Cl, O, TV)) dt. (44) 


According to the principle of extremal action we get Hamilton’s equation 
of motion for the true dynamical variables by varying dH with respect to 
(5 .|_, P- and their conjugate momenta mm . Also, we get one constraint 
equation from the variation with respect to the lapse, N. The equation of 
motion for the conformal factor is redundant. We do not consider it as a 
dynamical degree of freedom, but merely as a time dependent parameter. 

The constraint equation becomes 

pI+pI = ^ (« 2 + ^ (3) h) ■ (45) 


We are now in the position to quantize our theory. The canonical vari¬ 
ables are the two anisotropy parameters, /?+ and /3_, with their conjugate 
momenta and standard commutation relations. The conformal factor, fi, is 
treated as a c-number function that depends only on time. The constraint 
is imposed as an expectation-value equation, which in effect determines the 
time dependence of the conformal factor. The evolution of the quantum 
state is generated by the dynamical Hamiltonian, Hc> yn , which results in a 
Schrodinger equation for the wave function ip(/3+, /?_, t). 


JT I 

in— = H Dyn yj = 


N 


d 2 d 2 

+ 


+^ 

2m \d/3 + I 2N l 6 ) 


^ (46) 


The constraint equation produces a first order equation for 

^ = { pi +p l-^R) (47) 

We obtain a system of coupled integro-differential equations, (TIH1) and 
(HU), which exhibit an explicit dependence on the lapse, N. One requirement 
of any quantum theory of gravity is that measurable predictions should be 
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independent of the choice of lapse and shift functions, since different choices, 
after all, just represents different choices of basis. Shift functions have never 
appeared in this model, because the momentum constraints were satisfied 
identically, but lapse function is present. It is important to confirm that the 
system of equations given above give rise to dynamics which is independent of 
the choice N, despite its explicit appearance. To show this lapse invariance, 
let us note that, since the conformal factor is just a real function of coordinate 
time Q(t), one can find the inverse function (we assume that it is always can 
be done) t(Q), and use as the time parameter of the theory. In that case 
our wave function becomes i p(/3 + , f3-,Q). Furthermore, using 

_ d^dQ _ -dj)_ 

dt on dt on' 1 J 


we solve the constraint equation (1-T71) for fl 

N 2 ~ 

& = —(Co) (49) 

m z 

Here, Co = p\ + P- — R ■ In the last step we eliminate Cl from the 
equation for the wave function with the use of After simplification we 
obtain, 


— ih 


dip{0+,P-,n) 

dfl 


2\/(Co> 


8 2 d 2 

+ 


d(5 + d/3 _ 


+ m 




t/i(/3 + ,/3_,SJ). 


(50) 

This integro-differential equation is explicitly independent of the choice 
of the lapse function. In practice it is more convenient to deal with the 
system of the coupled equations where one can make any choice for N(t), for 
example, we take N(t) = 1. 

Another requirement of a reasonable candidate for a quantum theory is 
the existence of the classical limit. Basically that means that the theory 
should posses a regime which would mimic the results of the classical theory. 
The dynamics of the theory on the quantum level is governed by Huyn■ If is 
natural for one ask if the dynamics of the classical theory can be generated by 
the same Hamiltonian. Recall that the Hr) yn is defined prior the quantization 
by the standard Legendre transformation of the Lagrangian of the full theory 
with respect of the reduced set of dynamical variables. It is therefore logical 
to expect that if this Hamiltonian accomplishes the time translation on the 
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quantum level it will do so for the classical theory. It is not intuitively 
obvious that it is true for general relativity. It is certainly not the case for a 
general mechanical system, so the question is if there is such class of theories 
for which it is possible to reproduce the canonical equations of motion by 
defining the Hamiltonian over the reduced set of variables. 

In the case under consideration, the full canonical Hamiltonian, some¬ 
times called super-Hamiltonian, is Hq = NC° with the constraint given by 
(1381) . The canonical equations of motion are 


q a = {q a , H g }, and 

(51) 

p a = {p a ,H G }, 

(52) 

which subject to the constraint, 


C°(9„,p“) = 0. 

(53) 

Here, {q a } stands for {0, (3 + , /?_}. 

On the other hand, the equations of motion under the dynamical Hamil¬ 
tonian, H Dyn , are 

Qd = {qd,H Dyn }, and (54) 

P d = {p a ,H Dyn }, 

(55) 

and are subject to the constraint, 


C°(q d ,p d ,q e ,qe ) = 0. 

(56) 


In this case, {qd} = (/?+,/?_) and q e = Q It is shown (Appendix A) that the 
two sets of equations are equivalent. 

Actually we have proved the more general fact, that for a system of N de¬ 
grees of freedom described by a Hamiltonian of the type H = N a C a , where 
C a = 0 are the N c < N constraints of the first class in Dirac’s terminol¬ 
ogy, one can define the reduced Hamiltonian over N ( { — N — N c degrees of 
freedom. Moreover the canonical equations of motion generated by this re¬ 
duced Hamiltonian, together with the constraint equations, are equivalent to 
the full canonical theory. In turn, this means that, under these conditions, 
the constraints of the theory give rise to the equations of motion for the N c 
degrees of freedom. A detailed proof of this can be found in Appendix A. 
Here we just emphasize that the proof relies heavily on the fact that the 
constraints are of the first class. It seems almost trivial to generalize this 
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proof to the case of field theories of the same type, therefore the statement 
is true even in the case of full theory of general relativity. 

Now that we have established that the Hn yn works well as the Hamilto¬ 
nian of the classical theory, we just need to invoke Erenfest’s theorem. It 
is best formulated in the Heisenberg picture where time dependence of the 
operators are given by, 



Qd = H Dyn\, and 

(57) 


P^~{p d ,H Dyn }. 

(58) 

Based on the mapping between Poisson brackets and commutator established 
by {,} —^ ih[,] we see that, at least on the level of algebra, the Heisenberg 
equations of motion reproduce the classical theory. Therefore if we consider 
the states that are peaked around the classically allowed values of the canon¬ 
ical variables with small dispersion then 


• i * 

(Qd) — —J~{[QdiHi)y n ]), 

(59) 


(p) = ~([p d ,H Dyn ]), 

(60) 

together with 

( C°(q d ,p d ,q e ,qe )} = 0, 

(61) 


is guaranteed to give us the desired classical limit of the theory. 

As we have shown above, the states that can be called semi-classical are 
those that are highly peaked around a classically allowed configuration of 
canonical variables. Furthermore, such states will follow the classical trajec¬ 
tories in the phase space. However, are there new behaviors exhibited by 
those states that are far from being classical? The characteristic feature of 
the resulting quantum theory of the Bianchi IX cosmology is the constraint 
equation that is coupled to the Schrodinger equation. This equation gov¬ 
erns the evolution of conformal factor, Q(t). Although f }(t) is a classical 
variable with no dispersion allowed; nevertheless, it can deviate from the 
classical trajectory due to the coupling to the quantum state, iJj(/3 . + , /?_, t), of 
the system. For highly non-classical wave functions, one can expect to find 
substantial deviations or even qualitatively different behavior. It is known 
that the Bianchi IX cosmology does not predict a universe expanding with 
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an acceleration. Now in the light of this quantum theory we readdress this 
issue. 

It is easy to show that the quantum theory developed here is consistent 
with the deceleration behavior of the classical model in that it predict no 
accelerated expansion of the universe for any state of the system. This result 
is independent of the degree of non-classicity of the initial state function. To 
demonstrate this, we calculate explicitly the time derivative of the constraint. 


cl(C°) 

dt 


-^([C°,H Dyn }) + 


d(C°) 

dt 


0 


(62) 


The commutator of the constraint with H Dyn is zero. What we are left with 
is just the partial derivative with respect to time. It is convenient for our 
purposes to use the constraint equations in the form, 


= (63) 

m z 6 

Recalling the definition of the effective “mass” {m = 127rf?Qe 3n ), we find the 
time derivative of the constraint. 

N 2 

fi=-3— (pl+f_) (64) 

rn 


This explicitly shows that O < 0; therefore, the universe always decelerates. 


7 Conclusions 

The theory of canonical quantum gravity (Quantum Geometrodynamics) pre¬ 
sented here has been applied to the homogeneous Bianchi IX cosmology. It 
has been shown that the general formalism leads to a well defined and con¬ 
sistent quantum theory, which unlike other approaches outlined here (ADM 
and Wheelcr-DeWitt), allows one to make predictions relevant for cosmol¬ 
ogy. It has been proven that the equation of QG is independent of the choice 
of the time variable, (they are independent of the lapse function), and that 
in the classical limit, they agree with the classical Einstein’s equations of 
motion. This latter fact seems to be generalizable to the full held theo¬ 
retic case. The constraints are imposed as expectation-value equations and 
provide the coupling of quantum dynamical variables with the classical em¬ 
bedding parameters which, in turn, are responsible for the appearance of 
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the coherent classical spacetime. This way, the constraints of the classical 
theory have been given a new interpretation and physical meaning. On the 
practical side, the theory consists of a Schrodinger equation coupled to the 
first-order differential expectation valued equation. This system of equations 
can be successfully solved numerically, which was done in the case of vari¬ 
ous cosmologies including the Bianchi IX cosmology UU- It was originally 
expected that one might find effects different from the predictions of the 
classical theory for the highly non-classical states, for example accelerated 
expansion of the universe. At least for Bianchi IX model we have analytically 
established that his cosmological acceleration does not occur. The expansion 
of this model universe always decelerates. 
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9 Appendix A: A Proof of Equivalency of the 
Two Dynamical Pictures. 

For the sake of generality, consider a system with N degrees of freedom 
subjected to N c < N constraints of the first class with the full canonical 
Hamiltonian of the form 

H = N i C i (q a ,p a ), (65) 

where C l (q a ,p a ) = 0 are the constraints of the first class. Consequently, 
their Poisson brackets with each other weakly vanish, {C*,C'- 7 } = f l \C k . 
{q a }, a = 1 — N stands for the set of canonical variables, and N i} i = 1 N c 
are the Lagrange multipliers of the theory; the degrees of freedom that give 
us the primary constraints of the theory. 

Let us split the set of canonical variables, {g a }, cl = 1 — N, into two 
sets { q e }, e = 1 — N c and {gq}, d = (N c + 1 )tJV. We then define reduced 
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Hamiltonian by performing a Legendre transformation only over the second 
set of variables: 


Hr = p d q d (p d , Qd) - L(q d , q d (p d , q d ),q e , q e ) ( 66 ) 


Here the full Hamiltonian is, 


H = p i ild(p d , 'ld-P- qe)+p e <!e(p d , Pd-P- q e )~L(q d , ridip'pphp'd q e ), q, , q e (p d , q d ,p e , q e )). 

( 67 ) 

The pair of canonical equations of motion for the dynamical variables, 
{q d ,p d }, will be the same no matter what Hamiltonian is being used, simply 
because it is derived from the variational principle of the action, and this 
action is not altered by the Legendre transformation. Nevertheless, we show 
here in details that they are the same. 

The canonical equations of motion for {q d , p d } with the full Hamiltonian 
are, 

dH 

qd = W d ' and 68 

dH , . 

Pd = ( 69 ) 

Expressed in terms of the reduced Hamiltonian, they are, 


qd 


OHr 

dpd 


and 


( 70 ) 


. _ OHr 
Pd dq d ' 

What we need to show is that the right hand sides of these equations 
are the same. Equations m and (|70p just define the momenta. They are 
obviously the same because the momenta are defined with respect to the 
Lagrangian and not the Hamiltonian of the system. To see that the other 
equations are equivalent, we have to calculate the right hand side. 



dH R _ _ d dq d dL dhdq_ d _ dL_ _ dH 
dq d P dq d dq d dq d dq d dq d dq d 


( 72 ) 


Here we have used the definition of the momenta, p d = 

As far as the dynamical variables ({^d}) are concerned, it does not mat¬ 
ter which Hamiltonian is used. Incidentally, this leads to the vanishing of 
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the Poisson bracket of the full Hamiltonian, with the reduced Hamiltonian 
defined on the reduced phase space. It is important to note that the full 
Hamiltonian (H) is a constraint of the system, and therefore its Poisson 
brackets should vanish for consistency. 


fNMi 9 HdH R dll dii R 

R Dyn 


0 


(73) 


Let us shift our attention to the set of embedding variables {q e }- The 
canonical equations of motion for these are governed by the full Hamiltonian. 


• - dH - at d(Ji 
^ e dp e 1 dp e 


(74) 


Pe 


8H _ N ^ 
dq e 1 dq e 


(75) 


In the picture with the reduced Hamiltonian, we expect the constraint 
equations to provide us with the equivalent evolution. The constraint equa¬ 
tions are first order differential equations. One has to take the time derivative 
of the constraints to get the second order equations of motion. 


<7 

jC\q d ,p d , q eiP e ) = {&, H R } Dyn + — = 0 (76) 

We consider q e and p e , not as canonical variables, but as time-dependent 
functions, p e is related to q e through this definition. The time derivative of 
the constraint becomes, 


{C\H} Dyn + —f + —q £ = 0. 


dp e 


dq e 


(77) 


Therefore we have N c equations. We solve this system of equations for 
the N c time derivatives of the momenta. 


P k 


8C k \ 1 f9C > \ N dC l 
dpi ) \ dq e ) 1 dp e 



-l 

K". 


(78) 


Here, we have used the relation between the velocities and momenta ci e = 
ar l 

Note that the first term appears to be an incomplete Poisson bracket 
of the C l and C l over the embedding variables. We add and subtract the 
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appropriate term to complete these brackets, and we write down explicitly 
the content of the second term using the Hamiltonian H = N[C l . 


p k = - 


'dC k ' 

„ dpi ) 

( dC k ' 
V d Pi , 


-1 


-1 


'dc^ 
K dq e ) 
(dC l 


N, 


8C l (dC k 


-1 


N, 


dp e 

dC l 


dp e I dq, 


dPi J 

( 8C k ' 

l dpi 


N, 


'dC l dC l dC l dC l 


+ 


-i 


dq,j dp d dp d dqd 
'dC l \ AT dC l 

W ) Nl ~d<h 


Next, we combine the first three terms to produce the Poisson bracket of the 
constraints C l and C l taken over full set of variables {q e ,qd}- 


P k 



Ni{C\ C 1 } 


dC k \ 1 (d£\ dC l 
dpi ) \dp e ) 1 dq e 


(80) 


The Poisson brackets vanish in the weak sense, whereas the second term 
produces a Kroneker delta. 


p k « -8*N t 


dC l 

dq e 


= —Ni 


dC l 


dH 


dq k dq k 


(81) 


This result shows that, as long as constraint equations are satisfied for all 
times, they ensure the correct dynamical behavior of the embedding variables. 
The constraint equations can be used to generate the equations of motion 
for these variables which are equivalent to the canonical ones. 

To conclude this section, we would like to emphasize that the proof is 
based on the fact that the constraints of the theory are of the first class and 
the full canonical Hamiltonian of the theory is the constraint itself. These 
features are typical for general covariant systems. 


10 Appendix B: Invariant Basis for Bianchi 
IX Cosmology 


With the choice of coordinates we adopted, that of the three-sphere 0 < 6 < 
7r, 0 < p < 2i r, 0 < -0 < 4-7T, the Killing vectors are, 


6 = d vt 


(82) 


6 


cos pdg — cot 9 sin ipd v + 


sin p 

■ n 

sm 9 


and 


(83) 
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6s = - sin ( pd e 


.Q ~ . COS< ^0 

cot v cos g?oL H-— cL. 

sm 0 


(84) 


It is well known that the components of the metric tensor of a space- 
time which posses an isometry group take the simplest form in the so-called 
invariant basis, {e M }. The vectors of this basis are transported by Lie dif¬ 
ferentiation with respect to the Killing vectors. Mathematically, that means 
that 

[6,^ = 0 (85) 


To demonstrate the properties of the components of the metric tensor 
acquired in this basis, consider their Lie derivative with respect to the Killing 
vector 


&v) ( e AM T T 

= 9 {[£i,e y \,e v )+g{e, 1 ,[£i,e v ]) = 0 (86) 

In this basis, the components of the metric are constants along the curves 
generated by the Killing vectors of the given isometry group. In our case, 
there are three linearly independent Killing vectors, which means that the 
components of the metric do not depend on the three spatial coordinates but 
only on the fourth one - time. 

To generate the invariant basis for the Bianchi IX spacetime we proceed 
the following way. First, we choose a vector that corresponds to the time 
coordinate as the timclike vector of the invariant basis, eo = d t . This is 
possible because the Killing vectors, are the same for every hypersurface, 
£ t . Therefore, they are explicitly /-independent, and [d t , U] = [d t , = 

affdn = 0. The same time independence can be imposed on the spacelike 
vectors of the invariant basis by demanding that [eo, £j] = 0. To generate the 
rest of the basis we use, 


[6,e,] = [6,e™a m ]=0, (87) 

which results in a system of nine first order partial differential equations for 
three unknown functions e™ (9 , <p , ip) for each vector. One can solve these 
up to the choice of three free constants which corresponds to the freedom of 
choice of the basis vector at the “origin” (any point on the three-dimensional 
hypersurface). Using this freedom one can get the invariant basis that obeys 
the same commutation relations as the Killing vectors, namely [e,, ej\ = 
£ k ijek■ We list them here together with the dual basis of one-forms. 

e 0 = d t (88) 
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e.\ — — sin ibdg H- -t-d v — cot 6 cos ibd^ 

Sill 9 

(89) 

eo = cos + ——^-cL — sin ih cot 

sm 0 

(90) 

e 3 = dy 

(91) 

cr° = dt 

(92) 

a 1 = — sin i/jdd + sin 6 cos f)d^p 

(93) 

a 2 = cos ifd6 + sin 9 sin ij^dip 

(94) 

a 3 = cos 6dip + df) 

(95) 


If the timelike vector was chosen to be the unit normal to the hypersurface 
(£t), the metric would take the following form in this basis: 

g = —dt 2 + . (96) 
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